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Abstract 

We  consider  the  full  weak  convergence,  in  appropriate  function  spaces,  of  systems  of 
noninteracting  particles  undergoing  critical  branching  and  following  a  self  similar  spatial 
motion  with  stationary  increments.  The  limit  processes  are  measure  valued,  and  are  of 
the  super  and  historical  process  type.  In  the  case  in  which  the  underlying  motion  is  that 
of  a  fractional  Brownian  motion,  we  obtain  a  characterisation  of  the  limit  process  as  a 
kind  of  stochastic  integral  against  the  historical  process  of  a  Brownian  motion  defined  on 
the  full  real  line. 
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1.  INTRODUCTION 


In  this  section  we  shall  start  with  as  non-technical  as  possible,  and  therefore  occasion¬ 
ally  imprecise,  a  description  of  the  main  ideas  and  results  of  this  paper.  All  the  requisite 
qualifying  details  will  appear  in  the  following  sections. 

(a)  The  usual  setting.  Superprocesses  arise  as  infinite  density  limits  of  systems  of  branch¬ 
ing  Markov  processes.  We  start  with  a  parameter  ^  >  0  that  will  eventually  become  large,  a 
finite  measure  m  on  and  a  Poisson  point  process  H'*  on  with  intensity  measure  nm. 
The  K  =  0(/i)  random  points,  lo, . . . ,  ia"  of  the  Poisson  process  will  be  the  initial  positions 
of  a  system  of  particles.  Each  of  these  K  particles  follows  the  path  of  independent  copies  of 
a  Markov  process  Y,  until  time  t  =  If/i. 

At  time  l//z  each  particle,  independently  of  the  others,  dies  and  leaves  behind  itself 
a  random  number  of  “children”,  with  the  family  sizes  of  the  different  “parent”  particles 
independently  distributed.  The  mean  family  size  is  one.  The  individual  particles  in  the 
new  population  then  follow  independent  copies  of  Y,  starting  at  their  place  of  birth,  in  the 
interval  2//i),  and  the  pattern  of  alternating  critical  branching  and  spatial  spreading 
continues  until,  with  probability  one,  there  are  no  particles  left  alive.  It  is  clear  that  this  is 
a  purely  Markovian  system. 

The  basic  process  of  interest  is  the  measure  valued  Markov  process 

vti/ A\  _  {Number  of  particles  in  A  at  time /} 

A*  (A)  —  - - ,  (1.1) 

where  A  €.  8“^  =  Borel  sets  in  Note  that,  for  fixed  t  and  /i,  Xi  is  a  purely  atomic 
measure. 

It  is  now  well  known  that  under  very  mild  conditions  on  Y  the  sequence 
converges  weakly,  as  /x  — ►  oo,  and  on  an  appropriate  Skorohod  space,  to  a  measure  valued 
process  which  is  called  the  superprocess  for  Y.  If,  for  example,  V  is  a  Brownian  motion 
with  generator  A,  the  limit  process  is  the  super  Brownian  motion.  Details  can  be  found,  for 
example,  in  Ethier  and  Kurtz  (1986),  Walsh  (1986),  and  the  major  recent  review  by  Dawson 
(1993). 

(b)  A  non-Markovian  model.  The  motivation  for  this  paper  is  to  investigate  the  above 
problem  when  the  underlying  motion  of  the  particles  is  non-Mcirkovian.  This  problem  can 
be  divided  into  two  distinct  subproblems,  depending  on  what  happens  at  the  birth/death 
times  k/fi. 

The  simpler  of  the  two  cases  au-ises  when  a  general  non-Markov  motion  is  allowed,  but 
at  birth/death  times  each  of  the  new  particles  begins  a  motion  independent  of  those  of  its 
ancestors,  and  of  those  of  other  particles,  other  than  for  the  fact  that  it  begins  its  motion 
at  the  deathplace  of  its  immediate  parent.  Much  of  the  Markov  structure  is  thus  preserved. 
We  do  not  consider  this  case  in  the  present  paper. 

In  the  more  interesting  case,  the  one  which  we  do  consider,  the  particles  follow  the  talmu- 
dic  dictum  of  “Know  from  whence  you  have  come”,  (Talmud,  circa  400),  and  at  regeneration 
times  retain  some  memory  of  what  their  ancestors  did  in  the  past.  We  now  begin  with  a 
more  precise  setup. 
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We  append  to  the  state  space  a  cemetery  state  A,  and  adopt  the  convention  that 
^A)  s  0  for  all  functions  We  fix  0  <  p  <  oo,  and  recall  that  K  =  0(p)  is  the 

number  of  particles  alive  at  time  zero. 

In  order  to  label  our  particles,  define  the  family  of  multi-indices 

/  :=  {o  =  (Q!0)0!1  •••iQfjv):  oq  G  Af,  Oi  >  1,  »  >  1,  N  >  0}.  (1-2) 

Define  the  “length"  of  a  by  |a|  =  iV,  and  set  oj^  =  (oro, . . . , a,)  aind  a  —  i  =  (qq,  . . . ,  Q:|o|-i)- 
Induce  a  partial  order  on  /  by  setting 


-<  a 


0  =  q|,  for  some  t  <  |a|, 


and,  for  any  f  >  0,  write  q  ~  t,  if,  and  only  if 


M  <  <  <  (1.4) 

A* 

Furthermore,  for  two  indices  a  and  0,  set 

=  mjoc{A::  Oi  =  0i,  for  all  t  <  k},  (1.5) 

and  write  a  A  0  =  =  0\T<‘<fi-  Fo*'  ot  E  I,  let  {y’®(t),  t  >  0},  be  a  copy  of  the 

generic  process  Y  with  Vq  =  0,  chosen  so  that  the  family  satisfies 

y®(0  s  r®-‘(<),  (1.6) 

for  all  1  <  t  <  la].  Again,  for  each  a  €  /,  let  be  a  filtration.  We  require  that 

C  for  all  t  >  0,  whenever  a  <  0,  and  that  Y'*{t)  is  Tf  measurable.  In  general,  we 
think  of  as  the  o-algebra  generated  by  {y®(s)  s  <  t},  but,  even  in  the  most  interesting 
of  our  examples  to  follow,  this  will  not  be  the  case. 

Furthermore,  for  any  q,0£  /,  and  t  >  T®’^/p,  we  shall  require  that  Y°{t)  and  Y^{t) 
are  conditionally  independent  given  (When  ^  is  the  a-algebra  generated  by  V®, 

this  is  equivalent  to  requiring  conditional  independence  given  {Y°'^^{s):  0  <  s  <  T®’^/p}; 
i.e.  Children  recall  their  parent’s  behaviour,  but,  given  this,  ignore  their  peers.) 

It  is  not  hsu'd,  although  slightly  tedious,  to  show  via  a  standard  Kolmogorov-type  con¬ 
struction  as  in  Theorem  2.7.2  of  Ash  (1972)  that  such  a  system  is  well  defined.  This  almost 
defines  our  particle  system,  but  we  are  not  quite  finished,  for  we  have  yet  to  introduce  ran¬ 
dom  branching.  This  we  do  by  thinning  out  the  full  tree  of  particles.  For  each  q  €  /,  let 
be  an  independent  copy  of  a  generic  non-negative  random  variable  N,  which  represents  the  _ 
size  of  families  at  branching.  We  assume  throughout  that  EN  =  1.  Now  define  the  stopping 
times  ■  . 

fO  1  /»  cjf'c  tab 

^  J  mino<,<|«|  :  AT®!-  =  0|  if  this  set  0  and  ao  <  K  unannounced 

[  otherwise.  JusUf,cat;o’y 

3  By _ 

OtstribiJtion  ( 


DTIC  QUALITY  INSPECTED  S 


AviiiloD!  f'ly  Cotfi 

i  Avail  and /or 
'  Special 


We  am  now  finally  define  the  particle  paths  of  interest  to  us  by  setting,  for  each  q  6  /, 

V«(<)  =  |^ao  +  i“(0 


if  t  >  r®. 


Thus  Jf®  is  only  “alive”  (i.e.  ^  A)  at  time  i  if  it  has  a  continuous  stream  of  ancestors. 

The  definition  (1.1)  of  our  basic  measure  v^Jued  process  for  the  finite  system  is  now  more 
precisely  definable  as 

=  #W  G 

We  shall  also  find  it  useful  to  construct  another  measure-valued  process  associated  with 
the  above  branching  system  of  moving  particles.  Assume  that  the  generic  process  Y  has 
sample  paths  in  D[0,  oo),  the  space  of  S^-valued  cadlag  functions  (equipped  with  the  Sko- 
rohod  Ji  topology).  Furthermore,  for  any  function  / :  [0,oo)  and  t  >  0  define  a  new 

function  by  =  /(t  A  s).  Then  the  historical  process  for  the  above  particle  system 

is  well  defined  as 

HfiA)  = 

where  now  A  is  a  Borel  set  in  Z?([0,  oo),S^*). 

If  A  €  then  At  :=  {/  6  Z)[0,oo) :  f(t)  G  A}  is  a  Borel  set  in  Z?[0,oo).  Thus  it  is 
immediate  that  Xf  (A)  =  /ff(A«),  and  so,  in  principle,  knowledge  of  the  historical  process 
implies  knowledge  of  the  basic  process  Xf .  (The  converse  is  not  generally  true,  although  for 
special  processes,  in  high  dimensions,  it  is.  cf.  Barlow  and  Perkins,  (1993).) 

(c)  Self-similar  processes.  To  get  some  interesting  limit  theory  from  the  set-up  described 
above,  we  shall  have  to  make  some  assumptions  about  the  structure  of  the  generic  process  Y. 
In  particular,  we  shall  assume  that  {K(<),  f  >  0}  is  /f-self  similar  SJ**- valued  stochastic 
process  with  stationary  increments.  That  is,  for  any  n  >  1,  0  <  <  <2  <  •  •  •  <  h  >  0 

and  c>  0,  there  exists  a  real  H  such  that 


(y(cti ), . . . ,  Y{ct„))  k  (y(ti) . Y{tn)) 


(1.10) 


(/f-self  similarity)  and 


(y(<2  +  h)-Y{ti  +  h),  Y{t3  +  h)-  Y{i2  -b  h), . . . ,  Y{tn  +  h)-  Y{tn-i  -b  h)) 
=  -  Yih),Y{t3)  -  Y{t2),.. .  ,r(t„)  -  y(t„-i)) 


(1.11) 


(stationarity  of  the  increments).  The  common  alias  for  “/f-self  similar  with  stationary  in¬ 
crements”  is  “H-sssi”,  and  we  shall  use  it  in  the  sequel.  We  shall  also  assume  that 


0  <  /I  <  1,  and  J5||y(l)j|’'  <  oo  for  some  p  >  1/H. 


(1.12) 


It  is  worth  noting  that  mere  existence  of  the  first  moment  of  y(l)  forces  H  to  be  in  (0, 1],  and 
H  =  1  leads  to  a  degenerate  process  (Vervaat  1985).  Note  further  that  the  ff-sssi  property 
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of  Y  implies  that  £^11^(0  -  V'(s)||^  =  a|t  —  for  every  <,s  >  0  with  a  =  £||y'(l)||^,  and, 
of  course,  that  y’(O)  =  0. 

An  immediate  consequence  of  the  moment  assumption  (1.12)  is  that  Y  has  a  version 
with  continuous  sample  paths.  This  can  be  checked  in  a  variety  of  ways,  one  of  the  simplest 
being  via  a  version  of  Kolomogorov-Loeve  theorem,  (e.g.  BiUingsley  (1968),  Theorem  12.4). 
We  shall  assume  throughout  therefore  that  all  copies  of  Y  have  continuous  sample  paths. 
When  a  component  Yi  of  Y  €  is  a  one-dimensional  Gaussian  process,  it  is  the  well  known 
fractional  Brovmian  motion  and  so  can  be  represented  as 

Vi(0  =  t.  <  >  0  (1.13) 

for  H  ^  1/2,  where  ki  is  a  constant  and  B  is  a  standard  Brownian  motion  on  (—00,00). 
For  H  =  1/2,  (1.13)  is  not  really  very  informative,  and  Yi  is  simply  a  constant  multiple  of 
standard  Brownian  motion.  It  is  easy  to  check  (cf.  (1.15)  below)  that  the  covariance  function 
of  fractional  Brownian  motion  is  given  by 

Cov{Yiit),Y,{s))  =  ivar(Ki(l))[t2^  -f  -  \t  -  sp"].  (1.14) 

It  is  also  straightforward  to  check  that  every  (one  dimensional)  H-sssi  process  with  finite 
second  moments  has  the  same  covariance  function  (1.14)  as  the  fractional  Brownian  motion. 
However,  these  are  not  the  only  examples,  and  there  do  exist  non-Gaussian  one  dimensional 
/f-sssi  processes  with  finite  second  moments.  (In  fact,  with  all  moments  finite;  cf.  Surgailis 
(1981).) 

Our  main  weak  convergence  result  concerning  measure  valued  processes  (1.8)  is  stated 
and  proved  in  Section  2,  under  the  dependence  structure  outlined  in  Section  1(b)  above,  and 
the  assumptions  that  the  generic  particle  motion  Y  is  B-sssi,  satisfying,  together  with  the 
branching  mechanism  N,  an  appropriate  moment  condition.  The  result  parallels  the  weak 
convergence  to  the  superprocess  known  in  the  purely  Markov  case. 

(d)  More  on  fractional  Brownian  motion.  Since  the  general  convergence  results  just 
described  do  not  require  any  particulair  structure  for  the  puticle  motion  Y,  it  is  not  possible 
to  say  very  much  about  the  limiting  measure-valued  process.  To  do  so,  we  specialize  to  the 
particular  case  of  the  particles  moving  according  to  the  d-dimensional  fractional  Brownian 
motion,  (i.e.  each  component  of  Y  satisfies  (1.13)  for  some  ki,  and  the  same  H.)  Given 
the  general  terminology  of  superprocesses,  the  measure  valued  limit  process  could  well  be 
termed  a  “super  fractional  Brownian  motion”.  What  we  shall  establish  is  a  representation 
for  this  limiting  process  which  is  related  to  the  usual  super  Brownian  motion. 

To  do  this,  we  require  a  system  of  branching  Brownian  motions  on  (—00, 00),  so  as  to 
exploit  (1.13)  to  give  a  compact  and  useful  representation  of  the  particle  system. 

Retdning  the  indexing  notation  developed  above,  let  {{W^'*(0}(>o}a€/  ^  system  of 

branching  Brownian  motions,  all  stuting  at  the  origin.  Recall  that  at  t  =  0  there  are  K 
initial  points  for  this  system.  Let  W*,...,W^  be  K  extra  Brownian  motions  in  positive 
time,  independent  of  one  another  and  of  the  initial  system.  For  a  €  /  define  the  process 

B^it)  =  1[»>0](0W^“(0  +  l(»<0](0t^“'’(-0.  -00  <  t  <  00.  (1.15) 
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Then  each  is  a  Brownian  motion  on  the  entire  real  line,  centered  by  setting  B“(0)  =  0. 
Now  set,  for  each  a  6  / 

y“(t)  =  /‘  {(t  -  _  (_r)J-*/2jB“(dr),  t  >  0.  (1.16) 

00 

Finally,  define  the  thinned  system  {A^®}  as  at  (1.7).  However,  to  simplify  the  notation  for 
the  moment,  assume  that  all  the  initial  points  Xqq  &re  zero.  That  is,  all  the  paths 
start  from  the  origin  of  It  is  obvious  that  the  branching  system  so  constructed  fulfills 
all  the  requirements  of  the  previous  section,  with  the  filtrations  t  >  0  taken  to  be 
<r{jB*(s):  — oo  <  s  <  l}. 

Now,  however,  we  have  far  more  structure  than  we  had  before. 

To  see  this,  let  Hf  be  the  historical  process  for  the  system,  so  that  for  any  Borel  set 
y4ofC(»,»^)  and  t  >0 

=  -  E  MB”).  (117) 

Taking  the  integrated  form  of  the  measure  (1.8),  it  is  now  immediate  that,  for  any  nice  test 
function  — »  R, 

=  i  E  ■KX?) 

(It  is  easy  to  pass  from  line  to  line  here  if  one  keeps  in  mind  that,  with  the  exception  of  the 
one  stoch2istic  integral,  all  other  integrals  are  really  just  finite  sums.) 

In  Theorem  3.1  we  shall  show  that,  at  least  for  if  >  it  is  possible  to  pfiss  to  the 
fi  oo  limit  in  (1.18).  Since  the  final  line  in  (1.18)  is  then  a  functional  of  the  historical 
super  process  for  Brownian  motion  (albeit  on  all  of  R  rather  than  just  R+),  we  obtain  super 
fractional  Brownian  motion  as  a  function  of  super  (historic£d)  Brownian  motion,  giving  us 
the  main  result  of  Section  3,  and  title  of  the  paper. 

It  is  clear  that  the  explicit  representation  (1.18)  allows  for  a  far  clearer  understanding  of 
the  limit  process  than  does  the  more  general  procedure  which  we  descibed  earlier,  in  which 
nothing  other  than  the  existence  of  a  linut  process  is  established.  We  should  note,  however, 
that  the  previous  subsection  covers  a  much  wider  range  of  processes  than  mere  fractional 
Brownian  motion.  In  fact,  a  representation  like  (1.18)  would  seem  to  be  possible  only  for 
fractional  Brownian  motion,  or  other  processes  representable  as  functionals  of  systems  of 
branching  Markov  processes.  The  approach  taken  in  Section  2  works  in  far  greater  generality. 

As  we  shall  see  in  Section  3  below,  (1.18)  only  seems  to  work,  in  the  limit,  for  H  >  \. 
While  it  is  not  totally  clear  whether  this  is  an  artifact  of  the  proof  or  not,  we  believe  that 
the  result  may  not  be  true  for  smaller  H.  We  shall  explain  why  in  Section  3. 
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2.  WEAK  CONVERGENCE  FOR  GENERAL  SELF-SIMILAR  PROCESSES 


Throughout  this  section  we  shall  assume  the  branching  structure,  with  conditional  inde¬ 
pendence  of  new  particles,  described  in  Section  1(b).  We  shall  also  assume  that  the  particle 
paths  are  those  of  H-ssn  processes,  and  that  the  Poisson  distribution  of  initicd  points  is 
governed  by  a  finite  intensity  measure  m. 

Let  Aff(S^),  endowed  with  the  topology  of  weak  convergence,  be  the  space  of  finite 
Radon  measures  on  and  C([0yT),Mf{9t*))  and  Z)([0,T),Mf(Ji^))  the  spaces  of  con¬ 
tinuous  and  cadlag  functions  from  [0,  T]  to  When  T  =  oo  we  denote  these  more 

briefly  by  C(Mr(X^))  and  D(Mf(9t^)) 

The  following  is  the  main  convergence  result  of  this  section. 

Theorem  2.1  Suppose  that  for  some  p  >  2  both  £||V'(1)11*’  <  oo  and  EN^  <  oo.  If 
H  >  1/p,  then  the  sequence  of  Mp{fR^)-valued  processes  converges  weakly  in 

D{Mp{9i^)),  as  p.  —*  oo,  to  a  Mp{^i^)-valued  process  X.  Furthermore,  X  has  a  version 
with  all  sample  paths  in  C{Mp{ll^)). 

Remark:  We  shall  see  below  that  the  finite-dimension2d  distributions  of  A"**  converge  with¬ 
out  assuming  existence  of  any  moment  higher  than  second  for  the  branching  distribution 
and  for  every  H  G  (0, 1).  The  extra  assumptions  are  required  to  establish  tightness.  We  do 
not  know,  however,  whether  the  condition  H  >  1/p  is  necessary,  or  only  a  function  of  our 
method  of  proof. 

At  an  intuitive  level,  the  qualitative  aspects  of  this  requirement  are  readily  understood. 
Lack  of  smoothness  at  the  level  of  X*‘  can  be  due  to  one  of  two  sources;  irregularities  in  the 
particle  paths,  or  large  spurts  of  population  growth.  Recalling  that  the  sample  paths  of  the 
process  V  become  rougher  as  /f  decreases,  the  need  to  balance  the  moments  of  V  and  N 
against  the  smoothness  parameter  seems  more  than  reaisonable. 

Proof  of  Theorem  2.1.  The  proof  is  rather  long,  and  so  will  be  broken  up  into  a  number 
of  stages. 

Initially  we  shall  prove  that  the  sequence  converges  weakly  in  D(Mp(9t^))  to 

a  stochastic  process  A  which  has  a  version  with  all  sample  paths  in  C'(A//'(R‘^)),  where 
5?*  is  the  one-point  compactification  of  Only  at  the  very  end  shall  we  move  back  to 
D(Mp(9t^)). 

To  establish  the  weak  convergence  of  {A**},,  we  start  with  the  usual  reduction: 

Let  Cl  be  the  collection  of  read  valued  continuous  functions  on  which  have  a  continuous 
extension  to  the  one-point  compactification  of  R^  (that  is,  the  functions  with  a  finite  limit 
at  oo).  We  claim  that  it  is  enough  to  prove  that  for  every  tp  £  Ci  the  sequence  {A'‘(y>)}^ 
converges  weakly  in  I7(R)  to  a  process  which  has  a  version  with  all  sample  paths  in  C(R). 
It  will  then  follow  from  Theorem  3.7.1  of  Dawson  (1993)  that  the  sequence  {A^}^  is  tight  in 
D{Mp(9t^)).  To  establish  convergence  of  finite  dimensional  distributions,  note  that  since  the 
limiting  process  for  each  {A^(^)}^  has  continuous  sample  paths  and  R^  is  compact,  every 
limit  point  of  the  sequence  {A^}^  must  have  a  continuous  version.  It  then  follows  from 
the  assumed  convergence  and  Theorem  3.7.8(a)  of  Ethier  and  Kurtz  (1986)  that  for  every 
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i j  >  0, . . . ,  tp  >  0  and  €  Cj  the  sequence  ((<?), . . . ,  Aj^((^))}p  converges  weakly  in 
and  so  the  sequence  {(A/^, . . . ,  converges  weakly  in  (Aff (Si^))*’.  A  final  application 

of  Theorem  3.7.8(b)  of  Ethier  and  Kurtz  (1986)  shows  that  the  sequence  {A''}^  cannot  have 
more  than  one  limiting  point.  This  gives  the  full  weak  convergence  that  we  require. 

We  therefore  now  fix  6  C/  and  consider  the  sequence  {A**((^)}p  of  D(R)-valued  pro¬ 
cesses.  It  is  obviously  enough  to  consider  a  function  (p  with  a  bounded  first  derivative  (such 
functions  are  dense  in  C/;  see  Theorem  3.7.1.  of  Dawson(1993)),  and  henceforth  we  shall 
consider  only  such  functions.  We  start  by  establishing  the  convergence  of  finite  dimen¬ 
sional  distributions  of  this  sequence.  The  proof  is  a  reasonably  straightforward  application 
of  Lemma  3.1  of  Dynkin  (1991). 

We  shall  gi^e  enough  details  here  to  allow  the  reader  to  check  that  Dynkin’s  theorem 
applies  to  our  situation,  without  repeating  that  result  in  full  detail.  However,  it  is  worth 
noting  that  Dynkin’s  result  gives  convergence  of  finite  dimensional  distributions  for  br£inching 
systems  of  not-necessarily  Markovian  systems  in  an  extremely  general  setting.  His  result  is 
therefore  crucial  in  the  following  argument.  It  does  not,  however,  give  full  weak  convergence, 
(i.e.  in  the  function  space  setting  which  includes,  of  necessity,  tightness  of  some  kind). 
Furthermore,  in  the  non-Markov  setting,  it  works  only  at  the  level  of  the  historical  process. 

Choose  any  , . . . ,  tp  in  [0,  oo).  We  define  p  continuous  real  valued  functions  on  C(3?, 
by 

My)  =  <r»(y(<i)),  i  =  i,...,p. 

Observe  that  for  every  ft  >l 

(X, . «{(/,)) 

where  is  the  historical  particle  process  (1.9)  corresponding  to  A^.  Apply  now  Lemma 
3.1  of  Dynkin  (1991)  to  the  sequence  with,  in  his  notation,  K{dt)  =  dt,  branching 

mechanism  (p*^{x,z)  independent  of  t,fi  and  z  and,  finally,  with  the  additive  functional  A 
given  by 

A(5,a;)  =  X:i«.eB/.H 

i=l 

for  a  Borel  set  B  and  u  €  ^(R,??^).  The  application  is  straightforwud;  simply  treat  H'^ 
as  a  regular  (non-historical)  superprocess  with  respect  to  the  time  inhomogeneous  Markov 
process  Y  defined  by 

y(<)  =  (t,  y'(s),s<<). 

(A  detailed  example  of  how  to  do  this  is  given  in  the  following  section,  where  we  have  more 
need  for  rigour.)  Note  that  for  each  <  >  0,  Y{t)  is  x  C(3i,  8i**)-valued.  Since  Dynkin’s 
result  does  not  require  any  regularity  assumptions  on  the  Markov  process,  we  conclude 
immediately  that  the  sequence 

<W.(/l) . «{(/,)),  ,.  =  1,2,...} 

converges  weakly  in  W  as  p  —*  oo,  and,  therefore,  so  does  the  sequence 

{{X.'M,...,^'^),  ,.  =  1,2,...}. 
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This  concludes  our  treatment  of  the  finite  dimensional  distributions.  (Note,  by  the  way, 
that  we  have  not  been  able  to  say  anything  useful  about  their  structure!)  We  now  proce^ 
to  establish  tightness  for  the  sequence  D(X).  This,  in  itself,  is  a  rather  long 

calculation  that  needs  to  be  broken  up  into  a  number  of  steps. 

A  simple  extension  of  Theorem  7.2  of  Ethier  and  Kurtz  (1986)  (and  the  fact  that  is 
bounded)  show  that  it  is  enough  to  prove  that  for  every  e  >  0  and  T  >  0  there  is  a  6  >  0 
such  that 

(2.1) 

where  the  uniform  oscillation  function  u;  is  defined  by 

^(/,S,T)=  sup  1/(1) -/(s)|.  (2.2) 


As  a  first  step,  we  study  a  discrete  version  of  the  oscillation  (2.2),  which  is  somewhat 
easier  to  deal  with.  For  //  =  1,2,...,  S  >  0  and  M  >  1  consider 


u)^(f,S,M)=  .^max^  l/(il/M)  -  /(h/M)!- 


(2.3) 


We  shall  handle  this  oscillation  function  via  metric  entropy  methods,  for  which  we  need 
to  start  with  some  moment  calculations.  For  0  <  t'l  <  t2  <  A/p,  write 


o~h/M 

-  53  +  Xoo)l(t2/p  <  ’•(a)) 

-  53  v(j^“(*2/^)  +  Xoo)l(*2/^  <  ^(a))] 


:=  Ai(ii,t2)  +  A2(*i,i2), 


(2.4) 

where  for  a  ^  ti/p,  Aftt(ti,t2)  is  the  number  of  “descendants”  of  K®  alive  at  the  time  12/^. 
Observe  that 

A2(ti,t2)  =  53  +*ao)  -  V’(>'“(*2/a*)  +Xoo))i62/m  <  r(o)) 


N  {12 1  ft) 

==  E  Dj, 
>=1 


(2.5) 
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where  N{i2/ii)  is  the  total  number  of  particles  alive  at  time  t2/p. 

We  shall  use  two  simple  inequalities  which  caji  be  easily  checked  from  the  first  principles. 
The  first  is  a  version  of  Burkholder’s  inequality.  Let  Xi,X2,-  ■■  be  independent  zero  mean 
random  variables  such  that  E\Xi\^  <  oo  for  every  i  >  1.  Then,  letting  c  be  a  finite  positive 
constant,  in  general  dependent  on  p,  that  is  allowed  to  change  from  line  to  line, 

fii  >  m  iiy  m  V 

y=i  ,=i  j=i  ' 

for  every  m  >  1,  where  aj  :=  EXj.  Secondly,  let  Nn  be  a  critical  branching  process  starting 
with  a  single  individual,  with  a  progeny  distribution  having  a  finite  moment  for  some 
2  <  p  <  4  (say).  Then  there  is  a  c  €  (0,  oo)  such  that  for  any  n  >  1 

EN^  <  cnP-K  (2.7) 


Observe  that 


£|Aa(.,,.2)|'’<c(f:|  £  (D,-ED,)r+\EDtl'EN(i2/^)’’y  (2,8) 


Write 


VV{.2/m)  N{0) 


J=1 


k=:l 


where  5*  is  the  sum  of  all  the  terms  under  the  summation  in  the  left  hand  side  corresponding 
to  a  single  progenitor  existing  at  time  0  and  ./V(O)  is  the  initial  number  of  particles.  Then, 
by  (2.6),  we  have  that 


JVih/u) 


^^1  E  to  -  <  c(e(  E  Var(5,))"  +  £  E 

i=l  i=l  k=l  ' 

=  c(£;fV(0)'’/2(Var(Si))'’^^  +  £;7V(0)£;|5i|'’). 


(2.9) 


We  now  evaluate  the  moments  of  S\.  Let  be  the  size  of  the  tth  generation  of  the  jth 
particle  from  the  first  generation.  For  any  2  <  a  <  4  we  have 
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E\Si\*:=  E\Y,iDj  -  EDj)\‘ 
j=i 

<e(Ed\Z(Dj~EDj)\*J 

N}.  nL 


<  E  ■  ■  ■  E  -ED,,)...  (D,.  -  ED,,))) 

Si=i  i4=i 

.a/4 

<  ce( 53  - • .  1]  l»i/p  - 

Si*i  i4=i 


(2.10) 


=  cln/pL-i2/nr“E[N!,r 

<  c\ii/fi  - 

where  Eq  stands  for  expectation  taken  with  respect  to  the  Dj's.  Applying  (2.10)  first  with 
a  =  2  and  then  with  a  =  p  v/e  conclude  from  (2.9)  that 


7=1 


(2.11) 


To  complete  a  bound  for  (2.8)  we  need  to  evaluate  the  moments  of  N{i2lp)  appearing 
there.  Retaining  the  same  notation,  a  repeated  application  of  (2.6)  and  (2.7)  (in  the  notation 
of  (2.7))  yields  that 

/  MO)  ,  .  .\p/2  \ 

E{Nii2/py)  <  cE\Nioy  +  ( E  e{ni^  - 1)')  +  E  (2.^2) 

It  follows  now  from  (2.8),  (2.11)  and  (2.12)  that 

f^|A2(ti,*2)|'’  <  c\h/p  -  i2/p\^^.  (213) 

We  estimate  the  moment  of  A](ti,t2)  in  similar  manner: 

^iA,(n,t2)r<cp-^(  E  ^K-n-0)  +  g 

<  cp-''[(i2  -  iiY^^EN{iifpyf^  +  (*2  -  ti)'’-^£:A^(ii/p)] 

<  ~  *2//*!'’^*  +  1*1 //I  -  *2/^1'’“^)  > 

where  at  the  last  step  we  used  (2.12).  It  follows  from  (2.4),  (2.13)  and  (2.14)  that  for  any 
0  <  *1  <  *2  <  Mp 

E\Xf,i,W)  -  .»rf,,,(v>)l'  <  ^(lil//*  -  hb'f''  +  lil/c  -  ill  A’")-  (215) 
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We  can  now  use  our  moment  estimates  and  a  metric  entropy  argument  to  bound  the 
oscilkdion  function.  On  the  set  :=  {0, 1/p, . . . ,  A/  —  l/p,A/}  define  a  metric  d''  by 
setting 

For  an  e  >  0  let  N{t)  be  the  smallest  number  of  d'‘-balls  of  radii  not  exceeding  e  needed  to 
cover  the  whole  of  5**.  (For  ease  of  notation,  we  shall  not  explicitly  display  the  dependence 
of  N  on  p.)  It  follows  from  (2.15)  that 

In  particular,  this  implies  that  there  is  a  constant  K  <  oo  such  that  for  every  p  >  1 

<  K. 

It  now  follows  from  Theorem  11.6  of  Ledoux  and  Talagrand  (1991)  that  there  is  a  function 
h:  (0, 1)  — >  (0, 1)  such  that  h{S)  — >  0  as  d  — >  0  and  such  that  for  every  p  >  1  and  <5  G  (0,1) 

^  (216) 

An  application  of  Markov’s  inequality  now  establishes  a  discrete  version  of  (2.1):  for  every 
e  >  0  there  is  a  ^  >  0  such  that 

ns  p(a..(A-''(sp),«,T)  >£)<«.  (2.17) 

It  is  easy  now  to  bridge  the  gap  between  (2.17)  and  (2.1):  observe  that  for  any  t,s  £  [0,  T] 
we  can  place  t  in  some  interval  [ii/p,ti  +  1/p]  and  s  in  some  interval  [i2/Pi  *2  +  1//^]-  Then 
an  argument  very  similar  to  that  leading  to  (2.15)  shows  that  there  is  a  c  G  (0,  oo)  such  that 
for  every  /,  s  G  [0,  T] 


EWM  -  X'MV  <  c(|l  -  +  |l  -  »r),  (2,18) 

and  then  (2.1)  follows  as  before  by  an  application  of  Meu-kov  inequality.  This  completes 
the  proof  of  the  fact  that  {A’^}^converges  weakly  in  D[Mf{^))  to  a  stochastic  process  X 
which  has  a  version  with  all  siunple  paths  in 

It  remains  to  prove  that  the  limiting  process  X  is  supported  by  3^*,  in  the  sense  that  there 
is  an  event  fl+  of  probability  1  such  that  for  every  <  >  0  amd  u  €  f2+  we  have  Art({oo})  =  0. 

We  may  and  shall  restrict  t  to  be  in  the  intervad  [0,  T]  for  some  T  >  0.  For  a  fixed  r  >  0 
consider  a  bounded  Holder  function  on  3?*  such  that  0  <  <Pr(i)  <  1  for  any  x  G  3i^,  with 
'Pt{x)  =  1  if  ||x||  >  r,  (pr{x)  =  0  if  ||z||  <  r/2,  amd  with  a  Holder  constant  6,  say,  fixed. 
(5  >  0.)  (Such  a  function  always  exists  for  r  lau-ge  enough,  once  6  has  been  fixed.) 

Of  course,  (2.16)  and  (2.18)  apply  fully  to  tpr,  auid  since  the  constant  c  appearing  there 
depends  only  on  max^.^^  l¥’r(a:)|  and  on  the  Holder  constant  of  (pr»  we  conclude  that  for 
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every  f,s  €  {0,r],  and  any  r  large  enough,  (2.18)  holds,  and  the  constant  c  does  not  depend 
on  r.  Writing 

<('((,»)  =  (e\x;m  -  xiw,)r)''\ 

t,s  £  [0,7’],  we  estimate  next  the  d{f-diameter  Dr  of  [0,7’].  W'e  therefore  have 

<  c/*-' m>x,£|  E  l(l|i'“(l)  +  >  r/2)l(l  <  r(a))|>' 

<  52  *(l|V”(0ll  >  '•/4)l((  <  t(q))|'’ 


We  now  repeat  the  kind  of  arguments  that  led  to  (2.13)  and  (2.14)  to  see  that  for  any 
Q<t<T 

£|El(ll‘'“(')ll>'-/4)l((<T(a))|' 

S  c[£|  E[>(ll»''‘(')ll  >  ’■H)  -  f  (r“(i)ll  >  ’■/4))l(l  <  Ao))\” 

<  +  ^E\wy  +  i’(||r”(r)||  >  rnyEN(if^ , 

where 

IV  =  E  [l(r”'(<)ll  >  ’•/•f)  -  ^’(l|V»>(i)||  >  r/4)].  (2.19) 


Repeating  the  computation  in  (2.10)  we  easily  conclude  that  for  any  2  <  a  <  4 

E\W\‘  <  c,.-'p(||r«(0||  >  r/4)''‘  <  c;i->i>(||y"(r)||  >  r/4)“'’. 
Substituting  the  last  bound  into  (2.19)  and  using  (2.12)  we  conclude  that 

B\  E  l(ll>'“(<)ll  >  -•/4)l((  <  rCa))!'  <  6,i>'p(||y"(r)||  >  r/if'\  (2.20) 

a~t 

Furthermore,  let  A^r(O)  be  the  number  of  initial  particles  located  at  a  distance  of  more  than 
r  from  the  origin.  Since  Nr{0)  has  Poisson  distribution  with  mean  /im(R^),  where  Br  is  the 
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closed  ball  of  radius  r  centered  at  the  origin,  we  obtain  as  in  (2.12)  that  for  any  0  <  t  <  7" 


a~t  '  ' 

/  ,^r/4(0)  2\P/^  \ 

<c£(^Ar,/4(0)»’+(  E  E(N,,-1))  +  Y1 

A  consequence  of  (2.20)  and  (2.21)  is  that 

{dttMY  £  "['’(ll^"(J’)ll  >  ’■/*)"'*  +  ("■(«,'/,))''*]  :=  »(’■)• 

Observe  that  D(r)  |  0  as  r  -+  oo.  We  can  get  now  an  upper  bound  on  the  size  of 
supo<KT  l-^rCv^r)!  by  using  (2.18)  with  Theorem  11.1  of  Ledoux  and  Talagrand  (1992)  . 
We  have 

B  sup  W(s!',)-W(v.,)|<c /"‘'Vs- 
0«.«<T  Jo  \  J 


for  some  6  >  Q.  Therefore, 


=  cD{ry 


£  sup  |X,'‘(9,)|<£|.Vf(lS,)|  +  £  sup  W(^,)  - 

o<(<r  o^t.s^r 


<cifiENr^2(0}  +  D(rf) 

=  c(rniB^^2)  Dir^y 


(2.22) 


Note  that  the  constants  in  the  right  hand  side  of  (2.22)  do  not  depend  on  fi. 
Since  X^{(Pt)  =>  ^'(v’r)  weakly  in  Z?(3?)  as  ^  »  oo,  we  conclude  that 
sup  Xflipr)  =>  sup  Xt((pr)  in  distribution 

0<«<r  0<«<r 

as  fi  —*  oo,  and  so  by  Fatou’s  lemma  we  get 


E  sup  lX«(v9r)|  <c(m(5'/2)  +  £)(r)®). 
o<«r  \  / 


By  the  definition  of  ipr  we  immediately  obtain 


E  sup  A’«({oo})  <  c(m{B^j2)  +  D(r)*), 
o<«r  '  / 


and  since  the  above  holds  for  any  r  large  enough,  we  let  r  — »  oo  to  get 

E  sup  Xj({oo})  =  0, 
o<t<T 


implying  that  there  is  an  event  n+  of  probability  1  such  that  for  every  <  >  0  and  u  £  Q+ 
we  have  A’<({oo})  =  0. 

This  completes  the  proof  of  Theorem  2.1.  ■ 


3.  THE  FRACTIONAL  BROWNIAN  MOTION  CASE 

In  this  section  we  have  two  tasks.  The  first  is  to  carefully  define  the  process  that  we 
plan  to  use  in  obteuning  a  nice  representation  for  the  limit  of  branching  fractional  Brownian 
motions.  The  second  is  to  establish  the  weak  convergence  result. 


(a)  Historical  Brownian  motion  with  tails.  In  this  subsection  we  shall  closely  follow 
Perkin’s  (1992)  construction  of  historical  Brownian  motion,  and  of  certain  stochastic  integrals 
associated  with  it.  Much  of  what  we  have  to  say  could  well  be  swept  under  the  carpet  with 
a  totally  justihable  “it  is  easy  to  see  as  in  Perkin’s  paper”,  and  we  shall  actually  often  do 
so.  However,  we  nevertheless  require  enough  notation  to  properly  define  our  process,  and  to 
establish  some  technical  results,  in  the  following  subsections,  which  cannot  be  “easily  seen”, 
and  deserve  checking. 

Since  the  dimension  d  of  our  processes  is  fixed  throughout,  write,  respectively,  C  and 
for  the  space  of  continuous  functions  from  (— oo,  oo)  and  [0,  oo)  to  endowed  with  topology 
of  uniform  convergence  on  compact  intervals,  and  let  C  and  C"*"  denoting  the  corresponding 
a-algebras.  Let  ct  denote  the  canonical  filtration  on  C  (uncompleted  and  not  made  right 
continuous).  If  y  6  C,  to  6  C'*',  and  s  €  let 


(y/s/w)  =  I 


y{t) 

w{t  —  s) 


if  <  <  s, 
if  t  >  s. 


(3.1) 


Denote  Wiener  measure  starting  at  x  on  by  P*,  and  let  Mf[C)  denote  the  space 

of  finite  measures  on  C.  For  y  £  C,  set 


if  t  <  s, 
if  t  >  s, 


and  for  s  >  0  let 

MfiCy  =  {i/  €  Mf{C):  y  =  y^*l  for  v  a.a.  y}. 
Choose  V  €  M/-(C)(*)  -  {0}.  Define  €  Mf{C)  by 


PsA^)  =  (y/«/»")  e  A)  i/(dy). 

Let  ^  =  {(s,y)  G  R+  X  C :  y  =  y^*^},  and  for  (s,y)  G  ^  let  =  Ps,Sf  Finally,  define 
probability  measures  Ps,y  and  on  C  =  C([0,oo],^)  by 

KviA)  =  P..,({u;:  (s  +  (y/s/u,)(*+  ))  G  A}), 

P.AA)  = 
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where  v  6  — {0}  and  /I  €  C  =  Borel  sets  in  C.  If  Bt  (w)  =  ujt  denotes  the  coordinate 

map  on  C  2uid  {Cj+}  is  the  canonical  right-continuous  filtration  then  B  =  {C,C,Ct+,  Bt,  Pg^y) 
is  a  continuous  strong  Markov  process  with  a  Borel  semigroup. 

So  far,  we  have  done  nothing  other  than  to  repeat  Perkins’  (1992)  construction  of  the 
usual  historical  Brownian  motion,  with  a  redefinition  of  C  as  functions  on  the  whole  real  line 
rather  than  just  its  positive  half.  Further  details  can  therefore  be  found  in  Perkins’  paper, 
and  the  references  cited  therein. 

To  see  how  all  this  relates  to  the  constructions  of  the  previous  section,  we  require  some 
notation.  For  any  function  y  on  R,  define  y  by  y{t)  =  y{—t).  Simil^ly,  if  /I  is  a  Borel  set 
in  C(— oo,oo),  let  A  be  the  set  defined  by  .4  =  {y :  y  G  A}.  Finally,  if  m  is  a  probability 
measure  on  C(— oo,  oo),  define  m  by  m(A)  =  m(A).  We  leave  it  to  the  reader  to  check  the 
following  easy  lenruna. 

Lemma  3.1  In  the  construction  above,  take  v  =  Pq  at  (S.2);  i.e.  Wiener  measure  in 
"negative  time”.  If  Bt{Lj)  =  (t,Bt{u})),  then,  fort  >  0,  Bf  is  the  full  path,  from  — oo  to  t,  of 
the  “Brownian  motion  with  tails”  described  in  the  previous  section,  and  constructed  explicitly 
at  (1.15),  but  stopped  at  time  t.  We  let  P^p^  denote  the  law  of  Bt  on  C. 

Returning  to  the  general  setting  above,  it  now  follows  from  the  results  of  Fitzsimmons 
(1988)  or  Dawson  and  Perkins  (1991)  that  Bt  has  an  associated  Mf{^)  valued  superprocess 
H,  itself  a  continuous,  strong  Markov  process  with  a  Borel  semigroup.  If  Qo,^  denotes  the 
law  of  H  with  Hq  =:  Sq  x  v,  then  Ht~  Stx  Ht,  for  all  t  >  0  a.s.  We  let  denote  the  law 
of  H  on  where  fl/f  =  C([0,oo),  A!fjt'(C))  and  H  =  :  s  <  n).  The  process 

H  is  &  time-inhomogeneous  Borel  strong  Markov  process  with  Mp{C)  valued  paths.  This  is 
the  historical  process  that  we  shall  need  to  work  with.  In  the  setup  of  Lemma  3.1,  it  is  the 
“historical  Brownian  motion  with  tails”  of  the  previous  section. 

(b)  Some  properties  of  historical  Brownian  motion  with  tails.  Let  us  start  this 
section  by  agreeing,  henceforth,  to  denote  “Brownian  motion  with  tails”  by  BMT,  and  its 
historical  version  by  HBMT. 

We  now  need  to  collect  a  number  of  properties  of  HBMT.  The  first  is  a  straightfor¬ 
ward  application  of  Theorem  7.13  of  Dawson  amd  Perkins  (1991),  and  describes  the  weak 
convergence  of  a  particle  system  of  BMT’s  to  HBMT. 

Retadning  the  notation  of  Section  1(a),  let  be  a  system  of  branching  Markov 

processes  on  C,  with  transition  probabilities  as  at  (3.2).  Since  the  initial  vadues  B°'{0)  are 
all  of  the  form  (0,  B®),  where  B  is  a  BMT,  and  the  intensity  m  underlying  the  initiad  Poisson 
measure  H**  is  of  the  form  m  =  So  x  P^p^. 

Denote  the  associated,  particle  system,  historical  process  by  i.e.  for  A^C, 

Hf(A)  :=  /<-•  E  lA(^m  (3.3) 

a~t 

Note  that  fff  =  St  x  Hf,  where 

HtW  :=  /■-’  E  (3.4) 

Or~f 
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and  now  A  eC. 

By  Theorem  7.13  of  Dawson  and  Perkins  (1991)  it  follows  immediately  that  con- 
vergef  weakly  in  Z)([0,  oo),  to  the  sample  continuous  process  ff  defined  above.  An 

immediate  consequence  of  this  result  is 

Theorem  3.2  Under  the  above  conditions,  the  sequence  converges,  in  £)([0,  oo), 

Mf(C))  to  the  sample  continuous  process  H  defined  above. 


Our  next  step  is  to  show  that,  under  the  random  measures  Hf,  the  paths  y  behave  like 
the  branching  BMT’s  used  in  our  stochastic  integral  representation  of  fractional  Brownian 
motion.  This  is  crucial  to  all  that  follows. 

We  start  with  a  slight  extension  (because  of  tails)  of  a  somewhat  simplified  (because  of 
non-random  T  in  what  follows)  version  of  Perkins’  (1992)  Campbell  measures: 


Definition  3.3 
by 


Let  0  <  T  <  oo.  The  probability  measures  Pj-  on  {Clu  x.  C,H  xC)  defined 


Pr(A  X  B) 


E{\aHt{B))  ^  E{1aHt{B)} 

£{/fr(l)}  M 


(3.5) 


where  M  =  E{Ht{1)}  =  E{Ho{X)},  are  called  the  Campbell  measures  for  the  process  H. 


The  next  theorem  follows  directly  from  the  definition  of  Campbell  measures  and  the 
properties  of  //  as  in  the  proof  of  Theorem  2.6  of  Perkins  (1992). 

Theorem  3.4  Denote  a  generic  point  offljgxC  by  {u>,y).  Let  be  the  co-ordinate  mapping 
on  C.  Then,  for  each  T  >  0,  Bt{y)  :=  yt  is  a  BMT  stopped  at  T  on  {flu  xC,Hx  C,Pt)- 


The  importance  of  this  result  for  us  lies  in  the  following  theorem,  whose  proof  is  a  simpler 
version  (bar  the  extension  to  negative  time)  of  the  construction  of  the  Ito  integral  in  Section 
3  of  Perkins  (1992). 

Theorem  3.5  Fix  0  <  T  <  oo,  and  0  <  B  <  1.  Then  the  integrals 

con  be  defined,  Pj  a.s.,  in  the  same  fashion  as  an  Ito  integral.  Furthermore,  under  Pj, 
{Xi}  -oo<t<T  is  o  fractional  Brownian  motion  of  index  H. 

We  require  one  more  technical  construction  before  we  can  turn  to  the  main  result  of  this 
section. 

Definition  3.6  Let  0  <T  <  oo.  The  probability  measures  P^  ^  on  {flj{  xC  xC,HxC  xC) 
defined  by 

E{u  Ht(B)  IfT(C)}  E{1aSt(B)Bt(C)} 
Pt(AxExC).-  £{^,2(1)} - Mr  +  A/2 
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with  M  as  above,  are  called  the  product  Campbell  measures  for  the  process  H . 


Product  measures  of  this  form,  for  regular  super  Brownian  motion,  were  studied  in  Tribe 
and  Adler  (1993).  The  same  techniques  used  there  suffice  to  establish 


Theorem  3.7  Let  M  =  Denote  a  generic  point  of  (la  x  C  x  C 

by  iu,yi,y2).  Let  y,(f)  represent  co-ordinate  mappings,  and  fix  0  <  T  <  oo.  Then,  under 

Px  ,  the  pair  {Bi{t),  B2{t))  :=  iyi{t),y2it))  has  the  following  distribution: 


(B\(t)  B‘>(t))  —  i 


if  — oo  <  t  <  T, 
ift  >T, 


where  Wo  is  a  BMT,  Wi  and  W2  are  standard  Brownian  motions,  and  r  is  a  non-negative 
random  variable  taking  the  value  0  with  probability  Mf{T  +  M)  and  whose  conditional  dis¬ 
tribution,  given  {t  >  0},  is  uniform  on  [0,7].  All  four  variables,  end  r  are 

independent  of  one  another. 


(c)  Back  to  branching.  With  everything  properly  set  up,  w’e  can  now  return  to  the  setting 
of  (1.16)-(1.18),  where  the  branching  fractional  Brownian  motions  were  set  up  as  a  function 
of  a  system  of  branching  BMT’s. 

To  set  up  our  main  result,  let  be  a  branching  system  of  BMT’s  as  described 

above,  and  define  a  system  of  index-//  fractional  Brownian  motions  {V'®}o€/  ^  in  (1.16). 
As  pointed  out  in  Section  1(d),  this  system  satisfies  the  conditional  independence  structure 
we  require,  with  the  filtrations  generated  by  the  B°. 

Furthermore,  let  11^  be  a  Poisson  point  process  of  intensity  pm,  where  m  is  a  finite 
measure  on  5i^,  and  let  xi,X2,...  be  some  ordering  of  the  points  of  fl**.  Note  that  the 
a.s.  limit  of  is  simply  m.  Let  be  the  system  of  thinned  processes  based  in 

the  y®,  as  in  (1.7).  Finally,  let  A'**  be  the  particle  superprocess  (1.8),  and  let  be  the 
historical  process  (3.4)  corresponding  to  the  branching  system  of  BMT’s  described  above, 
with  weak  limit  H. 

Theorem  3.8  Let  be  the  particle  superprocess  described  above.  Assume  that  the  progeny 
distribution  has  p-th  moment,  for  some  p  >  2,  and  that  the  measure  m  has  all  its  mass 
concentrated  at  the  origin.  If  H  >  1/p,  then,  as  p  —*  00,  the  sequence  converges 

weakly,  in  D{Mp{^)),  to  a  Mp{lR^)-valued  process  X  with  a  version  with  all  sample  paths 
in  C{Mp{fR^)).  Furthermore,  X  has  the  following  representation,  in  law,  as  a  function  of 
the  historical  process  H  based  on  a  BMT: 

Xii'P)  =  J^v{x)Xt{dx) 

^  ~  -  {-r)+~^^^)y{dr))  Ht{dy), 

where  fit  ranges  over  all  bounded,  Lipshitz  functions.  The  inner  “stochastic  integral” 

is  that  defined  in  Theorem  S.5. 
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Remarks:  Note,  firstly,  that  if  H  =  1/2  the  fractional  Brownian  motions  are  regular 
Brownian  motions,  the  limit  process  is  super  Brownina  motion,  (3.8)  is  essentially  vacuous, 
and  so  there  is  nothing  new  to  study. 

Secondly,  Theorem  3.9  does  not  cover  the  case  when  the  initial  measure  is  not  concen¬ 
trated  at  the  origin.  In  this  setting  the  historical  process  H  has  to  be  changed  somewhat. 
In  particular,  rather  than  basing  the  historical  process  on  the  BMT  process  y,  we  need 
the  process  y{t)  :=  {z,y{t)),  where  z  is  a  random  variable  distributed  according  to  m  and 
independent  of  y.  In  that  case  the  representation  in  (3.8)  becomes 

~  (3.9) 

Note  that  for  all  Borel  A  €  SH,  Ht{C.,A)  =  m{A).  We  will  comment  on  where  the  proof 
needs  to  be  changed  to  accomodate  this  setting  when  appropriate. 

Finally,  given  the  structure  of  (3.8),  and  its  version  (3.10)  below  for  the  historical  process 
of  the  finite  system,  it  rather  looks  like  one  could  prove  Theorem  3.8  from  the  fact  that 
=»  along  with  a  simple  application  of  the  continuous  mapping  theorem.  However,  it 
follows  from  the  arguments  appearing  in  Nualart  and  Zakai  (1990)  that  the  mapping  from  the 
Wiener  space  of  BMT  realizations  to  defined  by  y  — ♦ 

can  be  extended  to  a  continuous  mapping  on  C  if,  and  only  if,  H  €  (1/2, 1].  Consequently, 
even  if  one  worked  to  get  this  approach  to  generate  a  proof,  it  would  work  only  for  /f  >  1  /2. 

Proof:  The  first  point  to  note  is  that  the  weak  convergence  and  continuity  of  the  limit  follow 
immediately  from  Theorem  2.1.  Thus,  all  we  really  need  establish  is  the  representation  (3.8). 
For  the  moment,  we  fix  t  >  0  and  bounded,  Lipshitz 

To  establish  (3.8),  we  commence  by  noting  that  (1.18)  gives  the  following  version  of  (3.8) 
for  the  finite  particle  system. 

L«^(x)xr(dx) 

r  .  (310) 

(Note  that  if  A''*(0)  ^  0,  as  assumed  here,  the  last  line  will  change  to  the  spirit  of  (3.9).) 

To  make  the  notation  somewhat  lighter  for  reminder  of  the  proof,  write  the  inner  integral 
in  (3.10)  as  F(y),  so  that  (3.10)  becomes 

XtM  =  J^v{F{y))  Hfm  (3.11) 

We  now  drum  that  there  exists  an  approximation,  say  F(,  to  F,  so  that,  for  each  fixed 
e  >  0,  and  each  ip  and  t  >  0, 

j^v{F({y))Ht{dy)  =►  J^fp{Ft{y))  Ht{dy)  as /x oo.  (3.12) 

Furthermore, 

M€^0,  (3.13) 
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and,  for  each  fixed  /i  >  0, 


j^v{F.(y))Ht(dy)ii  j^v(F{y))HHdy)  -  0.  (314) 

The  statement  of  the  Theorem  then  follows  via  £  standard  argument,  (e.g.  Billingsley,  1968, 
Theorem  4.2.) 

We  start  with  (3.12),  for  which  purpose  we  replace  the  innermost  integrand  in  (3.10)  by 
a  simple  function  from  (— oo,t]  — ♦  if,  of  the  form 

/(s)  =  if  s,_j  <s<sj,  j  = 

so  that  Ft{y)  becomes 

F.{y)  = 

y=i 

Define  a  map  M:  Mf{C)  — ♦  R  by 

>=i 

Note  that  this  map  is  continuous.  (Recall  that  C  has  the  topology  of  uniform  convergence 
on  compact  intervals).  Indeed,  the  function  C  — »  R  defined  by 

^(y)  =  v>(E  fiiyi^j)  -  yi^i-i)) 

>=1 

is,  obviously,  continuous  in  our  topology  and  bounded.  The  continuity  of  M  therefore  follows 
from  the  definition  of  the  vague  convergence.  Theorem  3.2  and  the  continuous  mapping 
theorem  now  suffice  to  establish  (3.12). 

For  the  proofs  of  (3.13)  and  (3.14),  we  require  some  notation  and  the  following  brief 
lemma. 

Let 

^2  =  if  ’l  <  <»}4 

^—00 

and  set  ||/||^,)  =  jLoo  f^i^)  Note>  for  later  reference,  that 

v,(r)  :=  (1  -  -  (r)?-''*  €  i|,  (3.15) 


for  all  t  <  oo,  if  if  €  (0, 1). 


Lemma  3.9  Let  tp  be  hounded  Lipshitz  and  f^g  €  L^-  Then  there  exists  a  universal  constant 
C  <  oo,  dependent  only  on  t,  such  that  for  every  ^  >  0 


/«»(*))  -  V>( /_%(>)»(*))]  (3.16) 
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whtrt 


IblUv  ==  ““P 


fb(j)-<r>(y)h 

^  ||x  -  y||  J 


is  the  usual  Lipshits  norm.  Furthermore,  (S.16)  abo  holds  when  H*‘  is  replaced  by  the 
limiting  hbtorical  process  H. 


Proof.  Note  that  the  left  hand  side  of  (3.16),  when  written  as  a  sum  rather  than  an  integral 
over  the  particle  historical  process,  is  equivalent  to 

cn^t  ^ 

X  [v>(  f_J(,)B>(d,))-V'[f_^9(>)B>(d,))\ 

<  ■£[#  of  pairs  (a, alive  at  time  (]^  •  (/(s)  —  p(5))B((is)J]^ 

<  ■  £|#  of  particles  Q  alive  at  time  tf  ■  ||v9|l*,p  •  £|  ^  I/(s)  -  p{s)l  B{ifa)|^ 

where  the  last  line  follows  from  (2.7). 

This  proves  (3.16)  for  the  case  of  finite  p.  In  the  case  when  //**  is  replaced  by  the  limiting 
historical  process  H,  the  left  hand  side  of  (3.16)  is  equivalent  to 

^[^{J  ^fis)w{ds))  -  ip[J  ^gis)w{ds))]Ht{dy)Ht{dw) 

<  lblli,>  -  9{s)]w{ds)])  Ht{dy)Ht{dw) 

=  ll^lli.>  -(Mt  +  M^)-  ( [  /_  J/(s)  -  g(s)]  y(ds)]  •  [ /_  J/(s)  -  .^(s))  u;(ds)] ) , 

where  Pp^  is  the  product  Campbell  measure  of  Definition  3.6,  and  M  =  £/fo(l)  =  m(5?^). 

Now  use  the  representation  of  the  pair  (y,  w)  given  in  Theorem  3.7,  to  evaluate  the  Pp^ 
expectation  above.  Incorporating  the  factor  (Aft  +  M^)  into  the  constant  completes  the 
proof  of  the  lemma.  ■ 

We  can  now  return  to  the  proofs  of  (3.13)  and  (3.14).  However,  these  are  now  quite  easy. 
Recall  that  in  order  to  establish  (3.12)  we  defined  by  approximating  t))e  function  (pt  of 
(3.15)  with  a  simple  function  on  S.  Call  this  approximation  ip\.  We  now  add  the  additional 
restraint  that  ||^(  —  — »  0  as  e  0.  In  view  of  (3.15),  this  is  easy  to  do.  Applying 

Lemma  3.9,  (3.13)  and  (3.14)  now  follow  immediately. 

This  completes  the  proof  for  fixed  t  and  fp.  To  complete  the  proof  in  general,  note  that 
we  can,  obviously,  repeat  the  above  argument  for  a  finite  number  of  times  and  test  functions. 
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and  so  show  that  the  joint  distribution  of  {A'/'j((^l), . . . ,  converges,  as  /i  — »  oo,  to 

the  joint  distribution  of  the  right  hand-side  of  (3.8)  with  the  same  <i, . . . ,  and  , . . . , 
Therefore,  (X/J, . . .  converges  weakly  in  Af/-(C")  as  /i  -*  oo,  and,  because  of  the 
tightness  guaranteed  by  Theorem  2.1,  Theorem  7.8  of  Ethier  and  Kurtz  (1986)  establishes 
the  full  weak  convergence  of  X**  the  right  hand  side  of  (3.8). 

This  completes  the  proof.  ■ 
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